function dx=Plant(t,x,flag,para)
dx=0;

iq=para(1);
ud=para(2);
w=para(3);

Rs=1.74;np=4;Ld=10;

dx(1)=-Rs*x(1)/Ld+np*w*iq+ud/Ld;